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ABSTRACT 


A summary .of modifications to Aerotherm’s Boundary. Layer Integral Matrix 
Procedure (BLIMP) code 'is presented c These modifications represent a preliminary 
effort to make BLIMP compatible with other JANNAF codes and to adjust the code 
for specific application to rocket nozzle flows. Results of the initial verifi- 
cation of the code for prediction of rocket nozzle type flows are discussed. 

For those cases in vjbich measured free stream flow conditions were used as input 
to the code, the boundary layer predictions and measurements are in excellent 
agreementc In two cases, with free stream flow conditions calculated by another 
JANWAP code (TDK) for use as input to BLIMP, the predictions and the datci V7ere 
in fair agreement for one case and in poor agreement for the other Cctse, The 
poor agreement is believed to result from failure of the turbulent model in BLIMP 
to account for laminarization of a turbulent flow. Recomiif^endations for further 
code modifications and improvements are also presented. 
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SECTION 1 


INTRODUCTION 

Boundary layer behavior along the walls of a rocket nozzle plays an im- 
portant role in the performance of the nozzle. The shear layer determines part 
of the thrust loss of the nozzle and the energy layer controls the heat trans- 
fer to the wall and the wall temperature. There is a well established need 
for a computer code v;hich can calculate boundary layer effects for flows v;ith 
large pressure gradients, chemical reactions, and a V7ide variety of v/all con*' 
ditions. To this end Aerotherm's Boundary Layer Integral Matrix Procedure 
(BLIMP) v;as selected by the JANNAP Boundary Layer Subcommittee as the standard 
performance evaluation method. 

The primary purpose of this report is to present the results of a pre-- 
liminary effort to provide a version of the BLIMP code that will serve as the 
standard boundary layer prediction tool for rocket nozzle flows. Section 2 con- 
tains a brief summary of the special modifications made to BLItiP to increase its 
utility to the rocket community, (Complete documentation of the original BLIMP 
code can be found in References 1 and 2.) The results of the verification tests 
for the four sets of data considered are contained in Section 3, In general the 
predictions are consistent v/ith the data in the wall j: eg ion of the boundary layer. 
This is based on wall heat flux comparisons, Thcnre are some discrepancies in the 
w^ake region; hsOv;ever, the predicted integral parameters are in reasonable agree- 
: ment with the data. In one case, for a laminarizing turbulent flovj, the hecit 
flux predictions did not agree with the measurements. Section 4 contains recom- 
mendations for further code developments and improvements. 



SECTION 2 


CODE MODIFICATIONS 

This section contains a brief discussion of the additions and modifica- 
tions to the BLIMP code which have been incorporated as a result of the present 
work* Details of these modifications are. contained in the interim user’s man- 
ual, Reference 3* This nev; version of the BLIMP code has been denoted as ver- 
sion J to distingvLlsh it from the standard ELIMP version^* Complete documenta- 
tion of the original code is presented in References 1 and 2* 

2.1 NAMELIST 

The BLIMP code was modified to accept the nozzle wall contour as (x, R) 
coordinates in the namelist format output by the Tv70- Dimensional Kinetic (TDK) 
computer program (Reference 4). These coordinates are then used to calcuJ.ate 
the nozzle wall length (assuming straight line segments between coordinate 
points) which is utilized in the internal coordinate system of the BLIMP code. 
The pressure ratio which is x^equired by BLIMP is also in the TDK neuiielist . A 
subset (max 40) of the input stations (max 500) is then selected by the user 
to be the BLIMP solution stations. 

2.2 THBRMOC HEMIC AL DATA 

The 35LIMP input format has been changed to accept the theriaocherriical 
constants for cal.culation of enthalpy, entropy, and specific heat In the poly- 
nomial form used by other JANNAF codes (Reference 5). 

2^3 INPUT/OUTPUT CHANGES 

In addition to the input changes described in Sections 2.1 and 2*2, the 
BLIMP code has been modified to accept input in either the International System 
of Units (SI) or the English Engineering System of Units* The entire output of 
the code has been formatted for SI units; however, in the case that English 
Units are used as input, they will also be used as output (the output headings 
remain in SI units). The outpvit has also been modified to include the boundary 
layer thrust loss (AF) , the total heat flux to wall (Q ) , and the total wall 
area calculated in the following manner: 
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where s is the wail length and q is the heat flux per unit area* 

An option has been added to calculate a new body contour or a nev7 inviscid. 
flov7 field contour based on a correction to the input contour* The new contovir 
is cjiven by 


^IKPUT ± '5* COS (j, 

where the + sign is for a new body contour and the - sign is for a new inviscid 
flov? contour* For example , if the input contour is the desired inviscid flow 
contour^ then the new contour , calculated with the H- sign, would be the cor- 
responding nozzle wall contour. The new contour can be punched on cards in a 
form suitable for input to the TDK program. 

2.4 REFIT OPTION 

. '/Ml option has been incorporated into the BLIMP code to cidju.st the nodal 
distribution as the solution moves dovjnstreara* This is accomplished by refit- 
ting the boundary layer profiles and shifting the nodes to insure that the nodes 
are always well placed* This option B>aintains a nodal distribution in the 
boundary layer which leads to. better accuracy in defining the velocity profile 
and more efficient use of computation time. 
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SECTION 3 


RESULTS OF WRIFIC/^TION ANALYSIS 


In a previous study (Reference 6) eKtensive comparisons of BLIMP predic- 
tions were made to experimental data. The cases considered included supersonic 
flow with zero pressure gradient, hypersonic flow v;ith zero pressure gradient 
and v;ith a favor^ible pressure gradient, hypersonic flow v;ith nonreactive blox^- 
ing, and superso 2 iic flow v;ith nonreactive blov/ing. Rather than repeat the broad 
range of cases considered in Reference 6, thei present study focused only on 
rocket nozzle type flows, i.e., those flows with very large favorable pressure 
gradients. One case considered in Reference 6 (Brott et al. , Reference 7) is 
of interest here and is discussed in Section 3.1. 

New data sources were sought which vjoiild provide detailed and accurate 
velocity and temperature measurements in the boundary layer, skin friction mca-^ 
surements, or heat flux measurements in flows with pressure gradients of the 
same magnitude as those in rocket nozzles. The pressure gradient similarity 
is desired because it is expected to be the dominant term in the boundary layer 
equations. In addition, the data should contain stagnation and free stream 
pressure measurements, and stagnation and v^all temperature measurements * The 
recent open, literature we^s examined and members of the JANNAF Boundary Layer 
Subcommittee were requested to provide any data sets of xvhich they had knowledge. 
Personal contacts at Rocketdyne and Aerojet v?ere also requested to supply use-- 
ful data on actual rocket firings. 

The Rocketdyne data discussed in' Section 3.2 is the only useful data 
which was found for actual liquid rocket engine tests. The input and some of 
the output for the predictions of this case are included in Reference 3 (the 
interim user’s manual for this version of BLIMP) as a sample problem. The hot 
air space shuttle main engine (SSME) model test data were also supplied by Rocket 
' dyne. These data consisted only of pitot tube measurements in the boundary layer 
at the exit of the nozzle; however, it was selected by NASA as a test case for 
^LIMP prediction. The most detailed, complete, and relevant data considered 
for verification was that of Back and Cuff el (Reference 8) and is discussed in 
Section 3.4. The input data and nozzle contours for the cases presented in 
Sections 3.2, 3.3, and 3.4 are contained in Appendix A. 



No adjustments to jthe BLIMP code were attempted to obtain agreement be- 
tween predictions and data. Modifications to the turbulent model and comparison 
with other models, e.g., Bushnell-Beckwith and Cebeci-Smith, should be investi- 
gated in future studies. Frequent reference is made to "the wall region" and 
"the wake region" of the boundary layer. In general these are rather loosely 
defined terms applied to turbulent boundary layers. The wall region refers ap- 
proximately to the 10-20 percent of the boundary layer near the wall. This is 
the region in which the "law of the wall" part of the turbulent model in BLIMPO’ 
is valid (see Reference 6 for a discussion of the turbulent model) . The wake 
region is the remainder of the boundary layer. 

3.1 DATA OF BROTT, YANTA, VOISINET, AND LEE 

The hypersonic flov; over a flat plate v^ith favorable pressure guradient 
data of Brottr et al. (Reference 7) is most nearly representative of 1102 ; zle 
flows of the data considered in Reference 6? although, the pressure gradient 
is not nearly as severe (more than an order of magnitude difference) . Repre- 
sentative comparisons from Reference 6 are shown in Figures 3-1 to 3-3 for stag- 
: nation conditions in air of 

== 1.013 X 10® N/m^ (10 atm) 

= 634. 5^K (610®R) 

The m omen turn thickness Reynolds nurrl;>cr (Re- “ p u G/u ) comparison in Fiqure 3-2 

u e e e 

provides information on the prediction of the overall boundary layer profile, 
since it is essentially an integral property c The skin friction coefficient 
^^f/2 ” comparison shov/n in Figure 3-3 Is a measure of how well the 

velocity profile is predicted near the wall. 

The conclusions from these three figures are that the predictions are 
very good at the wall (Figures 3-1 and 3-3) and that they are only slightly in 
error in the outer lav/ of the wall region (Figures 3-1 and 3-2). 

3.2 ROCKETDYNE: * 0^/H^ TWO-DIMENSIONAL NOZZLE 

This case* is representative of the type of data that can be expected 
^from hot fired nozzles using O^/H^ fuel systems. The gas side wall temperature 
and the wall heat flux distributions were measured; hov;ever, no boundary layer 


The data for this case were supplied by Hr. George. Osugi of Rocketdyne. 
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Figure 3-1. Linear-Log Velocity Profiled. Brott, et al . 

X/L - 5.583 (L = 0.3048m) [l = .0254m (1 in)] 
© Brott 's Data (u =•' 825 m/sec) 

BLIMP Predictions (u, = 823 m/sec) 










measureinents were made. Since heat flux is an important quantity and this case 
was for a representative liquid rocket fuel it was felt that predictions would 
be of interest. 

The' nozzle geometry, fuel mixture ratio, stagnation conditions, calcu-- 
lated axial pressure variation, and wall temperature variation v/ere provided, 
(The nozzle is .0762 m (3 in) wide.) The stagnation conditions were given as: 

- 4,6182 X 10® N/m^' (45.57 atm) 

= 3570^K (6430''R) 

M.R. - 6.15 


The injector pla 2 ie was at the entrance to the nozzle and the initially Iov7 heat 
transfer shown by the data (Figure 3--4) was assumed to result from the presence 
of a liquid layer near the injector. Accordingly, the prediction was started 
downstream of the injector and assumed to have an established boundary layer at 
the starting position. The mixture ratio in the boundary layer was assumed to 
he 6.15. The heat flux data and the results of the BI.IMP predictions are shown 
in Figure 3-4. The three cases shown in Figure 3-4 are for a fully turbulent 
boundary with a starting length of 0.0561 m, an initially laminar boundary layer 
with a starting length of 0.0561 m, and an initially laminar boundary layer with 
a starting length of 0.0286 m. (The starting length is the wall length from 
the start of the boundary layer growth to the first solution station. Adjust-- 
ment of this length has the effect of changing the boundary layer thickness at 
the first solution station. ) Transition for the Iciminar cases was at — 

400 and is sliovm on the figures. It is clear that the laminar boundary layers 
do not adequately predict the heat transfer. The turbulent predictions are ap- 
proximately 60 percent too large near the throat. Thus it seems that the bound- 
ary layer should not be either laminar or turbulent, and the possibility of a 
laminarizing boundary layer should be investigated. 

Schraub and Kline (Reference 9) have postulated that laminariza tion occurs 
in turbulent flov;s for values of the parameter K - exceeding 

about 3 X lO""® (see also Reference 10). The value of K at the first solution sta- 
tion is 1.33 X 10”**®. This strongly suggests that the turbulent boundary layer 
is. .laminarizing, which would result in a reduction in heat flux. The value 
of K increases through the tluroat (K^ == 2.81 x 10“^) indicating that the flow 
would become even more ** laminar like** in the throat region of the nozzle, thus 
reducing the heat flux from its turbulent value. 

This test case leads to two very important conclusions- First, it is 
necessary to have reasonable knov/ledge of the nature of the history of the boundary 
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Figure 3-4* Comparison of Measured and Predicted Heat Flux for 
Rocketdyne O^/H^ Nozzle Flow 

Throat radius: Rj = 3.93 x 10"^ m {0.155 in) 

Q*P (the test value of heat flux at the throat) = 

' 4.053 X 10^ watts/m^ (24.8 Btu/m"~sec) 

QPata 

BLIMP Predictions (S = 0.0551m, Transition Ro - 400) 

BLIMP Predictions (S^ = 0.0286m, Transition ReJ: - 400) 
BLIMP Predictions (S^ - 0.0561m, Fully TurbulelU) 
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layer at the first solution station (e.g*, did the boundary layer develop over a 
combusting liquid layer?) . Secondly, larainarization of the rocket nozzle bound- 
ary layer nay be a very significant effect and should be adequately modeled in 
the prediction tool. 

3.3 SSME MODEL TEST - EOT AIR 

These data from a heated air flow test in a 1:9.118 scale model of the 
space shuttle main engine wex“e provided by Rocketdyne* Pitot tube measurements 
of the boundary layer were made in the exit plane of the nozzle (area ratio = 
77.5). The pitot tube used in the wake region had dimensions 1.65 x 10”^ m 
O.D. (0.065 in) and 1.02 x m I.D. (0.040 in). Near the wall the probe was 

flattened with O.D. of 3.04 x 10’’^ m (0.D12 in). The relationships of these 
pirobe sizes to the bou nd a r y 1 ay er d i men i o ns are s how n in F i gu rr e 3 - 5 . ,S t ag n a - 
tion conditions V/ore given as: 


^ 1.027 X 10^ N/m^ (149 psia) 

‘ T = 634.5'‘K (610°R) 
o 

The wall contour v;as provided and, for the short duration of the test, the vjali 
temperature v^as assumed constant at 294.4'^K (530‘^R). The axial pressure distri- 
bution was provided from the output of the TDK program for the same conditions 
and geometry using the perfect gas option with j;atio of specific heats - 1.4 
and molecular weight - SB. SB. 

The results of the BLIMP predicfcio3is for boundary layer llach number dis- 
t r i b u t i o n at the ex it pi ane a x‘ e pr e s e n t ed w' i t h t li e data i n F i gur e 3 - 5 . T h e 
Mach numbers w^ere calculated from the measured pressures using normal shock re- 
lations. The major differences betv;een the predictions and the data are in the 
region 0.05 < y/6 < 0.37. Two BLIMP predictions were made, one for transition 
at Re^ -• 60 0 and the other for fully turbulent flow. There v;as no significant 
difference in the results at the exit plane for the two cases. The predictions 
for the transition at Re^ - 600 case are shovm in all the figures. No Icxminar- 
ization effects are expected since the value of K at the throat is 4,83 x 10'”'^. 

The predicted momentum thickness and displacement thickness are shown 
in 'Figure 3-6. The agreement with the values calculated from the data Is good; 
although, perhaps fortuitous. The predicted velocity profile and velocity pro- 
file calculated from the data are given in Figure 3-7, No temperature measure- 
ments were made in the boundary layer and all calculations from the data v/ere 


3-8 





1Q5, 6*/R 



Figure' 3-6. BLIMP Predictions for Momentum Thickness and Displacement Thickness 
SSME Model Test 

Throat radius: R_ = 0.014351 m (0.565 in) 

Momentum Thickness (G- test data result) 

Displacement Thickness (0: test data result) 






made by Rocketdyne assuming an adiabatic boundary layer.* Thus, calculated 
velocities would tend to be too large and this effect would increase as the 
wall is approached. This is reflected in Figure 3-7. The Mach number results 
are thus considered to be a more valid comparison. 

Based on the fact that tliese predictions were made from completely ana- 
lytical input, i.e., stagnation conditions and wall temperature, and pressure 
distribution calculated by the TDK program, the results are reasonably good. 

It should be pointed out that the BLIMP model for calculation of transport prop- 
erties (p and k) does not, in general, work well at very low temperatures such 
as those in the divergent region of this test case. In the present case, for 
example, the Prandtl number was calculated by the code to be 0.684 (one atmo- 
sphere, 273'’K) as compared to a more reasonable value of 0.73. The overall ef-- 
feet of these discrepancies appears to be slight, particularly since much of 
the boundary layer is turbulent. It is recommended, however, that a homogeneous 
gas option with special low temperature properties for air be added to the cooe 
to predict room temperature air flows. 

3.4 BACK AMD CUFFEL - SUPERSONIC NOZZLE WITH HEAT TRAl'iSFER 

These experiments were carried out in a cooled, conical nozzle with a 
convergent and a divergent half angle of 10° at the Jet Propulsion Laboratory , 
California Institute of Technology (Reference 8). The air flow was tripped 
well upstream of the converging section so that the flow was fully turbulent 
throughout the nozzle. These flow conditions are relevant to conditions in 
rocket engines operating at thrust levels for v?hich lariij-nanzatioii does not 
occur and provide a good test for the basic BLIMP capabd.3.ities ror irocket. nozz.le 
conditions. (The Urcoat value for K, defined in Section 3.2, is 1.0 x 10"®.) 

St.atic pressures r wall heat fluxes r and coolant-side wall temperatures 
were measured along the nozzle wall- Heat flux was determined by calorimetric 
measurements in circumferential coolant passages. Boundary layer surveys were 
made with a flattened pitot probe 0.000127 m (0.005 in) in heighth and with 
thermocouple probes. 

Velocity boundary layer profiles were measured at five stations; the 
locations of three of these are shovm in Figure 3-8. The edge values of for 
the predictions and the data are given in Table 3~1. In all the figures the 
test values of the various parameters asre used as normal.izing fa.ctors. The 
measured and predicted velocity profiles are shov/n in Figure 3-9 and the dis- 
p.lacement thickness (5*) and momentum thieJeness 0^ in Figure 3—10. The velocity 
profiles are in good agreement except for the position X/R^, = -6.414 v:here there 

^Personal communications v;ith Mr. Bill Wagner of Rocketdyne. 
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Figure 3-8. Schematic of Back and Cuffel Nozzle Showing Boundary Layer Probe Stations 

Rj = 0.0202 in (0.795 in) 
e = 9.89 

Nozzle Length (Converging-Diverging Sections) = 0.51 m (20.07 in) 




y/0 


Figure 3~9. -Measured and Pr^icted Velocity Boundary Layer Profiles For Back 
and Cuffel Data 

(6 values are calculated from data) 

(edge values of u are’ given in Table 3-'l) 

O Data 

BLIMP Predictions 




10^, 6*/R 



X/R^ 

Figure 3r10.. Momentum Thickness and Displacement Thickness For Back and Cuffel Data 
Throat radius: = 0.0202m (0.795 in) 

A - Data : 6* 

— BLIMP Prediction: 6* 


Data : 0 

■BLIMP Predictions:' 0 



are differences (less than 10 percent^ however) in the wake region. The inte- 
gral properties 6 and 6* which reflect profile shape are in good agreement ex- 
cept at the position X/R^ = 10, 797, The difference in is a consequence of 
a difference in the predicted and measured temperature profiles. The total 
temperature^ profile is shov/n in Figure 3-11, (The normalizing factor is the 
value of T* - T- from the data,) Since the velocity profiles at X/'Rm ^ 10.797 
are in very good agreement , the differences in total temperature can be attrib- 
uted to differences in static temperature. These slight differences ma.y result 
from the use of a constant value of 0,9 for the turbulent Prandtl number. The 
different shape of the static temperature profiles will cause different density 
profiles. The density-velocity product, mass flux, shown in Figure 3-12, is 
used for the calculation of 6*, It is the large negative contribution to 6* 
from the region above = 1 the data that causes the difference in 

the measured and predicted values of 6*, 

The near wctll behavior can be evaluated by examining the heat f.lux data 
shov 7 n in Figure 3-13, The excellent agreement supports the conclusion that the 
near wall region is properly accounted for in BLIMP, The more detailed prose- 
tation of the velocity profile at X/R^ - 10,797 shown in Figure 3-14 shows that 
the data does not extend into the viscous sublayer region; therefore, one must 
rely on the heat transfer measurements for compciring wall behavior. The excel- 
lent overall agjreement in velocity px-ofiles is apparent 

TABLE 3-1 

EDGE VALUES OF VELOCITY FOR BACK A® CUFFEL DATA 



-11.21 

-6.414 

10.797 

i!^(m/s). Data 

42.0 

79.4 

1100.2 

u^(m/s), BLIMP 

42.4 

80.0 

1111.6 


3,5 SUMMARY 

The two cases of RocKetdyne data represent almost completely analytical 
predictions, i,e,, from given stagnation conditions and wall geometry an axial 
! pressure varition was predicted and used in the boundary layer predictions. 

The two-dimensional nozzle predictions do not fit the data at all; hov^ever , 

1 

■;’tlie"‘f low could well be of the laminarizing turbulent type for which the turbu*-' 
Aent model in BLIMP J has not been verified. Thus there is a need to investigate 


"Total temperature, T", is the measured stagnation temperature profile in the 
boundary layer. 
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Figure 3-11. Measured and Predicted Total Temperature Profiles far 


Back and Cuffel Bata, X/R^ = 10.797 
(0 calculated from data = 0.652 x 10"‘*m) 

Q Data (T* - T* = 476“K, T^ = 35S°K) 

-- — BLIMP Predictions (normalized by 476°K, T* - T* = 492'’K, T = 356°}(1 

^ , vv . . w ^ 
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the modeling of laminari nation in the code. The SSME air flow data are reason- 
ably vjell predicted. More detailed data of this type would be useful for code 
verification. In contrasty the data of Back (and of Brott) contained measured 
pressure distributions v;hich were used as input to BLIMP, Tlie resulting predic- 
tions shov; much better agreement with the data. 
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SECTION 4 


RECOMMENDATIONS 

As a result of working with BLIMP and rocket nozzle problems the follow- 
ing additions and modifications to the present version are suggested, 

1. In order to provide the capability of comparison of results and 
potentially for the selection of the "better” model, the eddy 
viscosity mociels of Cebeci-SmitJi and Bushne Hr: Beckwith should be 
added to the code, Fi^rthex^ verification of the present model (of 
Kendall) and comparison of the nev7 models using experimental flows 
providing tests of various features of the models which are important 

i to nozzle conditions should be made* Of particular interest is the 
-■prediction of laminar ization, 

2, Cold flow model tests using air, or other homogenous gases, could be 
treated by the BLIMP code without the use of the complete chemistry 
solution. ' Simplified procedures and improved computation economy 
could be achieved through the addition of a homogeneous gas option. 

This option could also offer the advantage of improved lov; temperature 
thermodynamic properties a.nd transport properties. 

. 3. For greater compatibility v;ith the T^'? 0 “Dimensional Kinetic program 
the follow'ing recommendations are made: 

e. Addition of Kinetics 

© Improved matching of edge conditions on velocity, pressure, 
and their first derivatives, and improved curve fits of 
pressure g3:adient* The present method vjas not designed for 
the type of pressure distribution encountered in nozzle flov/ 
problems. For rocket nozzle flovrs the largest terms of the 
boundary layer momentum equation contain the gradient of 
pressure and the gradient of edge velocity ;; thus accuracy in 
these terms is most important. 

4. The complexity of the input can be decreased by an expanded use of 

namelist input and built-in default values. This would greatly aid the * 
inexperienced user and pose no liinxtations to the generality of the code 
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5. Overall accuracy and computer usage could be improved by the follow-- 
ing: 

© . Expanded three temperature range to provide for better lov; 
temperature curve fits of thermodynamic properties. 

© Binary diffusion option for economy for problems with complex 
chemistry, 

& Improved transport property calculations, especially for 
light elements and air. 

o Improved built--in first guess of fully turbulent flows. The 
built-in first guess is for a laminar p^rofile. For those 
problems v/ith v;e 11- established turbulent profiles at the first 
solution station, a turbulent first guess v/ould improve coraputer 
economy. 

o Improved convergence criteria tailored to noz 2 le flows. 

The present convergence criteria have been developed and tested 
for reentry conditions v/here pressure gradients are not nearly 
so severe as in rocket nozzles. 
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APPENDIX A 


INPUT DATA AND NOZZLE CONTOURS 



A.l ROCKETDYNE 0^/H^ TWO-DIMENSIONAL NOZZLE INPUT DATA 
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Figure A-l, Bodtetdj/ne 2"l) flozzle Contour 

ArG0 Ratio *■ 4,0 

Aspect Ratio (Width/R,,) •» 22.8 

Length = 0,1484 m (5.B46 In) 

R-f *> 0.003937 m (0J55 In) 
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2,029445+00, 
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2,614028+00, 
2,7429071-00, 
2,876083+00, 
3,013656+00, 
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3.302420+00, 
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3,037030+00, 
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1,152272+00, 
1.252649+00, 
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1,4 66244+00, 
1.560619+00, 

1 ,6557811-00, 

1.753664+00, 

1.854484+00, 

1.950426+00, 

2.065478+ 00, 

2,176077+00, 
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2.406418+00, 
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5,710910+00, 

6,309415+00, 

6.9891831-00, 
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1.071532-100, 

1 . 168681+01', 

) .269742 + 00, 
1.375037+00, 

1 .483150 + 00, 
1.576250+00, 
1.671950+00, 

1 .770269 + 00, 
1.871546+00, 

1 .9759431 00 , 
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2,1948/5+00, 
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3 . 377480100 , 

5.531009+00, 
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4.000961+00, 
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4.806527+00, 

5 . 281012 + 00 , 
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1.000000+00, 1, 000047100, 
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1,019912+00, 1,022879+00, 1,026139+00, 
1,042648+00, 1,047810+00, 1,053470+00, 
1,032501+00, 1.091040+00, 1,099451+00, 

1,133157+00, 1.141589+00, 1,150005+00, 

1,133745+00, 1,192209+00, 1,200684+00, 
1 ,234f.77+00, 1,243220 + 00, 1, 251750 + 00, 
1,256120+00, 1.294720+00, 1.303574+00, 
1,336174+00, 1.346905+00. 1,555690+00, 
1,391032+00, 1.399916+00, 1,400856+00, 

1,444040+00, 1.453695+00, 1.463010+00, 

1 ,593831 + 00, I i.5141 46 + 00, 1 .524521+00, 
1.566319+00, 1,576906+00, 1,587531+00, 

1,630467+00, 1.641379+00, 1,652239+00, 
1,696427+00, 1.707643+00, 1,718835+00, 
1,764364+00, 1.775910+00, 1,757460+00, 
1,834436+00, 1.646335+00, 1,858272+00, 

1,906799+00, 1.919074+00, 1,931427+00, 
1,978109+00, 1.980467+09, 1,998645+00, 
2.040762+00, 2.0S135H 00, 2,061931+00, 
2.104617+00, 2.115391+00, 2,126224+00, 
2,169814+00, 2,130746+00, 2,191781+00, 
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2,3/3714+00, 2.385446+09, 2,397253+00, 
2,444873+00, 2,456739+00, 2,466842+00, 
2,517524+00, 2.529863+00, 2,542113+00, 
2,591906+00, 2.604501+00, 2,617145+00, 
2,668009+00, 2.660866+00, 2,693307+00, 
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2,825439+00, 2,638847+00, 2,852287+00, 
2,906672+00, 2,920441+00, 2,934110+00, 
2,939599+00, 5.003635+00( 3,017649+00, 
3,074196+00, 5,036484+00, 3.102834+00, 
3,160413+00, 3.174955+00, 3.189650+00, 
3,246262+00, 3.263031+00, 3,277916+00, 
3,337769+00, 3.368207+00, 3,399030+00, 
3.527437+00, 3,560757+00, 3.594396+00, 
3,734159+00, 3.7/0306+00, 3.007013+00, 
3,958049+00, 3,997024+00, 4,036284+00, 
4,198639+00, .4,240380 + 00, 4.26262 0 + 00, 
4,457331+00, 4,502545+00, 4,548471+00, 
4,739769+00, 4.789453+00, 4,839613+00, 
5,046424+00, 5.099019+00, 5,151910+00, 
5,366631+00, 5.421082+00, 5.475932+00, 
5,700368+00, 5.757903+00, 5,815929+00, 
6,053776+00, 6,114419+00, 6,175748+00, 
, 6,424513+00, 6.487705+00, 6.551^60+00, 

i 6.807819+00, 6,872217+00, 6,937106+00, 
i 7,196123+00, 7.260897+00, 7,325375+00, 
7,581593+00, 7,644993+00, 7,70/951+00, 
7,955182+00, 8.015377+00, 8,074924+00, 
0,305030+00, 8,360216+00, 8,414355+00, 
8,616/37+00, 8,666398+00, 6,712597+00, 
PnA8(l) = l ,0, 

,9993289+00, .9968/73+00, 

,9560030+00, ,9299963+00, ^8910000+00, 
.7699944+00, ,6719799+00, .4690063+00, 
.3629402+00, .3478866+00, ,3330965+00, 


1,632/16+01, 1.664208+01, 1.696096+03, 
1.839654+01, 1.866032+01, 1.902024+01, 

2,053477+01, 2.093249+01, 2.133799+03, 
2.304184+01, 2,348920+01, 2,352260+01, 
,725721, 1,627636, 1,460524, 1.354854, 

1,012336, 1.00057, 

1.000193+00, 1.000446+00, 1.000814+00, 

1,005579+00, 1.004635+00, 1,005354+00, 
1,012576+00, 1.014779+00, 1 , 0 1 72 1 6 + 0 0 , 
1,029758+00, 1.033641+00, J . 037939 + 00 , 

1.059689+00, 1.066530+00, 3,074078+00, 
1.107339+00, 1.116284+00, 1,124712+OC, 
1,158446+00, 1,106860+00, 3.175316+00, 
1,209154+00, 1.217668+00, J. 226151+00, 
1,260545+00, 1.268688+00, 1.277465+00, 

1,312056+00, 1,320719+00, 1.329436+00, 
1.364493+00, 1.373298+00, 1.382158+06, 
1,417809+00, 1.426776+00, 1.455601+00/ 
1,473170 + 00, 1. 465347 + 00, 1 . 493607 + 0 6 , 
1.534S78+00, 1.545325+00, 1.555S21+0C, 
1,593175+00, 1,608909+00, 1.619669+00, 
1 .6632.'i9 + 00r 1.674271 + 00, 1,685293+00, 
1,730139+00, 1.74 1516+00, 1.75288/1 + 00, 
1,799123+00, 1,610843+00, 1.822583+00, 
1,870521+00, 1,882406+00, 1.894547+00, 
1.943911+00, 1.956545+00, 1.967740+00, 
2,009293+00, 2,019763+00, 2,030227+00, 
2.072535+00, 2.083214+00, 2.093959+00, 
2 . r3.7030t00, 2.147900 + 00, 2,158833+00, 
2,202379+00, 2,213971+00, 2.225091+00, 
2,?.70078+00, 2,201507 + 00, 2,292804 + 00, 
2.333777+00, 2.350390+00, 2. 362116+00, 
2,«090'+9+00, 2.420913 + 00. 2.432842+00, 
2,480972+00, 2.493124+00, 2.505266+00, 
2. 554497+00, 2. 5f, 6963 + 00, 2,579419 + 00, 
2.629742+00, 2,642468+00, 2.655301+Ou, 
2,706748+00, 2,719/04+00, 2.732/78+0O, 
2.783495+00, 2,798679+00, 2.612014+00, 
2,8oS3«3+00, 2.879432+00, 2.892986+00, 
2,947925+00, H.96S873+00, 2.975655+00, 
3.031693+00, 3.045855+00, 3,060038+00, 
3,117105+00, 3.131513+00, 3,146051+00, 
3.204135+00, 3,218/96+00, 3,233553+OU, 
3.292/84+00, 3^307685+00, 3,322713+00, 
3, 430466+00, 3,4'62328 + OOr 3,49/(592+00, 
3.625521 + OOr 3 .6<, 3422 + 00, 3.696580 + 00, 
3,8/14056 + 00, 3,001509 + 00, 3.919508 + 00, 
4.076077+00, 4,116400+00, 4,157254+00, 
/(, 325402+00, 4,3667(30 + 00, 4,/4l2766 + 00, 
4.595142 + 00, 4,642591 + 00, /(, 690797 + 00 , 
4,890759+00, 4 ,9422 1 } + 00 , 4,994121+00, 
5.203168 + 00, 5,258715 + 00, 5. 5.12512 + 00, 
5,531205+00, 5.586984+00, 5,643586+00, 
5,874363+00, 5.933782+00, 51993590+00, 
6.237332+00, 6.299329+00, 6,361835+00, 
6.614366+00, 6,679095+00, 6.743151+00, 
7,001699+00, 7,066552+00, 7.131360+00, 
7.389/72+OOf. 7,454065+00, 7,517084+00, 
7,7/0439+00, 7,832556+00, /, 894100+00, 
6,133751+00, 6,191821+00, 8,248870+00, 
0,467353+00, 0,519138+00, 8,569623+00, 

3,757247+00, 6.600224+00, 6,5033o7+C0 

,9883995+00, ,98.13572 + 00, 

,8344955*00, ,0029922+00, 

,3939462+00, ,376279/1 + 00, 

.3185643+00, ,3042776+00, 


A- 4 



130. 

131. 

132, 

133, 
33^1, 
133. 
136, 

1 3 V , 

136., 
13<). 
1«(0.. 
HU 
142 , 

143., 
144, 
14 5. 
146. 

147., 
146. 


14 V. 

150.. 

151. 

152. 

153. 

154. 

155. 

156. 

157. 
156. 

159, 

160 . 
Ifcl, 
162. 

163. 

164. 

165. 

166. 
167, 

166., 


169, 

170, . 

171, 

172, 

173, 

5.74., 

175. 

176., 
3 77. 

176. 
179, 
100 , 
10 ). 
162. 
163, 
104. 

165, 

166. 
167, 
166, 
169, 

190, 

191, 

192, 

193, 

194, 

195, 

196, 


,29022VOtOO, 

,2234505+00, 

. 1626601 + 00 , 

, 1067770+00, 

,6313536-01, 

,4976724-01 , 

,4845610'-01 , 

.4722082-01, 

.4569290-01, 

,4462136-01, 

,4336092-01 , 

.4209341-01, 

.4056591-01, 

,3930432-01, 

,3769626-01, 

,3646531-01, 

,3478407-01, 

,3339100-01, 

,3202026-01, 

,30672t>0-0] , 

,2909004-01, 

,2779675-01, 

.2653674-01, 

.2329023-01, 

.2434995-01, 

.2319425-01, 

,2231226-01, 

,2121962-01, 

,2039652-01, 

.1937512-01, 

,1336933-01, 

.1765257-01, 

.1673634-01, 

,1605959-01, 

,1521266-01, 

,1459472-01, 

,1361535-0 I f 

, 1325422-01, 

,1254018-01, 

, 12053.33-0 1, 

.1138176-01, 

. 1 092617-01, 

.1049259-01, 

,.9924930-02, 

,9535710-02, 

,6637295-02, 

.7910104-02, 

,7,086801 -02, 

,6435813-02, 

,5326943-02, 

,5355159-02, 

,4869621-02, 

,4428248-02, 

,3993379-02, 

,3661789-02, 

,3309336-02, 

.3063052-02, 

.2790575-02, 

,2615975-02, 

,2364562-02, 

,2207600-02, 

,2071597-02, 

. 1929964-02, 

.1763566-02, 

.1660161-02, 

.1547700-02, 

SEND 


,2764128+00, 

.2107971 + 0.0, 
,1513146+00, 
,9593494-01, 
,5516107-01, 
.4947361-01, 
.4640110-01, 
,4707443-01, 
.4578011-01, 
,4452357-01, 
.4300399-01 , 
.4173674-01, 
.4046910-01 , 
.3694769-01, 
.3750164-01 , 
,3607424-01 , 
,3466451-01 , 
.3337528-01, 
.3164240-01 , 
,3030111-01, 
.2696557-01 , 
,2744120-01 , 
,2621366-01, 
.2524974-01 , 
.2406606-01, 
,2291699-01 , 
.2204437-0) , 
.2095966-01 , 
,2014424-01 , 
,1913061-01, 
, 1336696-01, 
, 1742405-01, 
.1672241-01, 
,1584703-01, 
,1520462-01, 
.1439767-01, 
.1301313-01, 
.1307236-01, 
,1254319-01, 
,1186567-01, 
,1136942-01, 

, I077260--01 , 
.1034460-01, 
;978a270-02, 
,9264265“02, 
.85J7966-02, 
,7662924-02, 
.6977633-02, 
.6334363-02, 
,5731663-02, 
,3262225-02, 
.4796453-02, 
.4356923-02, 
,5916921-02, 
.3605672-02, 
,3306768-02, 
.3011104-02, 
.2760556-02, 
.2546929-02, 
,2370067-02, 
,2191073-02, 
,2009000-02, 
.1869390-02, 
.1761031-02, 
.1644464-02, 
.,1522649-02, 


,2628242+00, 

,1983850+00, 

,1402307+00, 

,8943256-01, 

,5070962-01, 

,4913162-01, 

,4605607-01 , 

,4672511'>01 , 

,4542807-01 r 

,4416952-01 , 

,4290566-01 , 

,4137928-01, 

.4011234-01, 

,3861503-01 , 

,3737501-01, 

.3566336-01 , 

.5427786-01 , 

,3289362-01 , 

,3153108-01 , 

.2993072-01, 

,2862208-01., 

,2734140-0) f 

,2591793-01.- 

.2496065-01, 

,2378537-01, 

,2283647-Oi, 

,2177793-01, 

,2093349-01, 

,1989354-01, 

,1910794-01, 

.1813416-01, 

. 1719763-01 , 

,1650346-01, 

.1563653-01, 

,1500148-01, 

.1420247-01 f 

.1362502-01, 

,1289213-01, 

. U36980-01 , 

,1170025-01, 

. 1 123012-01, 

. 107B267-0)., 

, 1019865-01 f 

,9797675-02, 

,91.43167-02, 

.8260701— 02, 

,7548961-02, 

,6869901-02, 

,6235711-02, 

.5636894-02, 

.5169365-02, 

^4703356-02, 

.4251115-02, 

.3837061-02, 

.3530771-02, 

,32536)7-02, 

,2939805-02, 

,2710681-02, 

,2533571-02, 

,2304422-02, 

.2127079-02, 

,1990027-02, 

.1648887-02, 

.1738138-02, 

,1620599-02, 

.1516203-02, 


,2494665+00, 

, 1662244+00, 
,1294374+00, 
.»026920"0l , 
,5035799-01, 
,4376817-01, 
,4791657-01, 
,4659512-01, 
,4532550-01, 
,4381375-01, 
,4254953-01, 
,4128026-0 1 , 
.3975515-01, 
,3341962-01, 
,3695226-01, 
,3556165-01, 
,3415992-01, 
,3251237-01, 
,3115596-01, 
,2932396-01, 
,2526004-01 , 
,2698666-01 f 
.2567528-01, 
,2467323-01, 
,2374997-01, 
<2261593-01, 
,2151346-01, 
,2067 6 5 0-01 , 
,1964445-01, 
, 1386674-01 , 
,1790107-01, 
,1716185-01, 
,1628619-0 1 , 
.1562673-01, 
,1479933-01, 
.1419567-01, 
.1343864-01, 
.1289567-01, 
,1219318-01, 
. 1170664-01, 
„ i 107261-01 f 
. i0t,3l 1 3 ”01, 
,1005458-01, 
.9658920-0?, 
,9021 952-02, 
.8143439-02, 
,7435734-02, 
,6762945-02, 
,6133379-02, 
,5542659-02, 
.5076380-02, 
,4611)705-02, 
,4165120-02, 
.3758561-02, 
,3456608-02, 
,3160939-02, 
.2932064-02, 
.2698626-02, 
,2465431-02, 
,2280339-02, 
,2109'429-02, 
,1970550-02, 
.1827309-02, 
.1715018-02, 
. 159654(.-02, 
,0000000 , 


,2363410+00, 

,1743224+00, 

.1109491+00, 

,7151720-01 f 

,5010706-01/ 

,4679959-01, 

.4756982-01 , 

,4624500-01 , 

,449/600-01 p 

.4371584-01 , 

,4219238-01 , 

,4092337-01, 

<3966070-01, 

,3829036-01, 

,3658935-01, 

,3517260-,01 , 

,3377541-01 r 

,3239881-01 , 

,3078171-01 , 

,2945657-01 , 

.2615797-01 , 

.2688823-01 , 

,2558200-01 r 

.2438749-01, 

,2347142-01 , 

,2234306-01 , 

,2148132-01, 

.2042105-01, 

.1962106-01, 

,1662703-01, 

,1788262-01, 

.1695820-01, 

,1607093-01, 

,1541882-01, 

, 1460015-01, 

. 1400612-01 , 

,1325410-01, 

.1271603-01, 

.1220253-01, 

. 1154331-01, 

,1100152-01, 

. 1048144-01 , 

, 1006729-01 , 

.9521635-02, 

. 8756996-02, 

,3026902-02, 

,73226’93-02, 

,6538962-02, 

.6034801-02, 

.5446765-02, 

,4983142-02, 

.4518904-02, 

.4080842-02, 

.3ottl569-02, 

,3302521-02, 

,3088505-02, 

,2800762-02, 

,2628761-02, 

.2451545-02, 

.2272856-02, 

.2090706-02, 

.1950369-02, 

. 1806022-02, 

. 1691699-02, 

,1572334-02, 

,0000000 , 


i 
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190.. 
.199, 
200^. 
2 9'U 
20^, 

20/U 
205^ . 
206. 

207, 

208, 

209, 

210 , 
211 , 
212, 
213, 
21^, 

215, 

216, 

217, 

218, 

219, 

220 . 
221, 
222e. 

223. 

229, 

?2S. 

226. 

227. 

226. 

229, 

230, 
251. 

232, 

233. 

239.. 

235. 

236. 

237. 

258. 

259. , 
290. 
29 U 
292. 


12 


0.0 ■ , 

002 

,006 

,01 

,025 

,06 

.« , 

7 

1.0 

1.5 2 

.5 


10 .V5 

1 

3 0.5 




0.0 

.05 

.12 

,25 

,35 

,45 

.75 

,010351 

,85 

.95 

,98 

1.0 


1,02732+06 

0,15650+04 






0.0 






.04 11 

,023 ,013- 

,9 

,9 

600, 

3.3900 + 05- 

.z 

N NITiiOGEN 


14,000 -,766 




0 oxygen 


16,0 -.232 




100, 


1000, 5000 

t 



N2 C0!.0 Alf? 


NONEi N 2 


. G 

iOO, 5000 

3,51515 





-1054.545 


6.06115' 

3', 51515 
-1054,545' 


6,08115 

02 COLO AIR 


none; (7 2 

G 

100, 5000 

3,53 







-1059.09 

6,06ae9 

3,53 




-1059,09 

6^06a89 


,15 

• 6 


299,9 

299,9 

29zj^£| 

29^, a 


29^1,9 
29a, ^ 
299./? 

29a, a 


29a. a 
29a. /j 
29a. a 
. 29a, a 


29a. a 
29/4, a- 
294.a 
29a. a 


29a. a 
29a. a. 
29a. a 
29a. a 


2?a.a' “ 
29a. a 
29a. a 
29a. a 


2 9a. a 
29/4. a 
29a. a 




oaioo A 

oaaoi a 

0a?02 A 

oa^oo A 
oaaol A 
oaao2 A 
05100 A 
07100 A 
07200 A 
07300 A 
08100 A 
09600 A 
nioo A 

11201 A 

11202 A 
13100 A 

I 


13LAST ' 

15201 A 

15202 A 

15203 A 
1520a A 

16201 A 

16202 A 

16203 A 
1620a A 
16601 A 
16602, A 
16603 A 
166oa A 

16605 A 

16606 A 

16607 A 

16608 A 

16609 A 
3 66.10 A 
1 6 6 3 A 
16612 A 

A 
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X: w ru w ru 




10 15 20 25 

X/R^ 


SSME Mode;! Nozzle Contour 
= 77.5 

376 m {14,825 in) 

51 ni (0.565 in) 



A. 3 


INPUT DATA FOR BACK AND CUFFEL PREDICTIONS 


1, 

10100^^202 

1021200000G 

JPL DATA 

AIR CAROS 



01100 

2, 

3,,. 

4. 

2 

2/J 

0*01 

0*05 

0-1 

0,2 

0*9 

l.O 

2,0 

02100 

03100 

03201 

5» 

•93,761 

3 . a 1 7 

6*107 

4.234 

9 * 36 1 

0*^80 

4,575 

03202 

6, 

/|*609 

i^.633 

i|,64l7 

4,678 

9*720 

0*781 

-4,90 7 

03203 

7. 

e, 

o; 

5*063. 

13 

0*0 

5*250 

0*002. 

5*379 

0.006 

0,01 

0*025: 

0*06 

0,15 

03204 

04100 

04201 

10,. 

oj\ 

0*7 

UO 

1,5 

2.5 

4,0 


04202 

11,. 

12, 

10 0,95 
0,0 

1 3 .5 

0*05 

0*12 

0,25 

0,35 

0*45 

0,6 

04300 

04^101 

13,. 

0,75 

0*B75 

0«95 

0.985 

0,99 

1*0 


0 4 4 02 

14. 

15, 

1.0 

*20033 

*20833 

*20833: 

,20033 

*20333 

,20633 

,20033 

05100 

05201 

16, 

*20633 

• 1925 

.1592 

,1374 

,H53 

*09325 

,07833 

05202 

17, 

*07233 

• 06617 

*06665 

,06700 

*0745. 

*06517 

,10708 

05203 

18, 

19. 

20, 
21, 
22, 

23, 

24, 

25, 

26, 
27, 

.134 ,1685 ,169- 

. 10,21 
240,9?- 

0,0 . . 

0,44 11,02 0,018 

36,74 
2 

R MITROGEN 14,000 -,760 

D DXYGEN 16.0 -,232 

100, 1000, 5000, 

C,9 

0t9 



05200 

07100 

07200 

07300 

08100 

09600 

lUOl 

11201 

11202 

13100 

20*. 

N2 COLD 

AIR none; N 2 


G IDO* 

sooo. 




29.. 

30, 

31, 

32, 

33 , 
3'n 
35 , 
30, 
37 , 
30, 

39, 

' 10 , 

01.. 
0 2 , 
03 „, 
00, 

45 .. 
46, 

47 .. 

40 .. 

49 .. 

50 . 

51 , 

53, - 

54, 

55 .. 

56 , 

57 , 

50 .. 
59 , 

60 .. 

01,. 


3,51515 

»L054,545 

D2 AIR 

3,53. 
-1059,09 


6,00115 
NCINE: 0 


6,06459 


3,51515 

»1054.545 


3,53 

-1059,09 


6,00115 
C 100. 


6,06409 


5000, 


J 

J 

J 

J 

J 

J 

J 

J 

J 

J 

J 

J 

J 

J 

J 

J 

J 

J 

J 

J 

J 

J 

J 

J 

J 

J 

J 

1 

2 

3 

4 


13UAST 


0.9972 

0.9972 

0*9972 

0*9972 

0,9972 

0*9972 

0*9972 

iSiOl 

J 

0.9972 

0.99835 

0^99257 

0*98703 

0 , 97 <l 0 

0-93704 

0,6659 

15102 

J 

0,7807 

0,02767 

0,6760 

0,016555 

0*4987 
0«0121 1 

0 * 35 1 2 

0,2317 

0^1489 

0^05927 

15102 

15201 

15202 

15203 
. 15204 

J 

J 

J 

J 

0 

690, 

690, 

690^ 

690 R . 

690, 

690o 

690* 

16201 
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